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Abstract 

We consider the particle spectrum and event shapes in large N gauge theories in 
different regimes of the short-distance 't Hooft coupling, A. The mesons in the small 
A limit should have a Regge spectrum in order to agree with perturbation theory, 
while generically the large A theories with gravity duals produce spectra reminiscent of 
KK modes. We argue that these KK-like states are qualitatively different from QCD 
modes: they are deeply bound states which are sensitive to short distance interactions 
rather than the flux tube-like states expected in asymptotically free, confining gauge 
theories. In addition, we also find that the characteristic event shapes for the large A 
theories with gravity duals are close to spherical, very different from QCD-like (small 
A, small N) and Nambu-Goto-like (small A, large N) theories which have jets. This 
observation is in agreement with the conjecture of Strassler on event shapes in large 
't Hooft coupling theories, which was recently proved by Hofman and Maldacena for 
the conformal case. This conclusion does not change even when considering soft-wall 
backgrounds for the gravity dual. The picture that emerges is the following: theories 
with small and large A are qualitatively different, while theories with small and large N 
are qualitatively similar. Thus it seems that it is the relative smallness of the 't Hooft 
coupling in QCD that prevents a reliable AdS/QCD correspondence from emerging, 
and that reproducing characteristic QCD-like behavior will require genuine stringy 
dynamics to be incorporated into any putative dual theory. 



1 Introduction 



Confining gauge theories play a central role in particle physics: they provide the real the- 
ory of strong interactions (QCD), and possible models of dynamical electroweak symmetry 
breaking (EWSB) and dynamical supersymmetry breaking. With the help of the anti-de 
Sitter/conformal field theory (AdS/CFT) correspondence many of the qualitative features 
of confining gauge theories can be explained, and warped extra dimensional theories have 
been used to parametrize strongly coupled, almost conformal, models of dynamical EWSB. 
The AdS/CFT analogy [1] has been pushed so far as to suggest that most properties of 
QCD can be understood even quantitatively using a simple Randall- Sundrum (RS)-type ex- 
tra dimensional model called AdS/QCD: a warped extra dimension with UV and IR cut-off 
branes, and bulk fields appropriate for the chiral symmetries of QCD [2]. While many results 
show a surprisingly good agreement (0(5 — 10)%) with the QCD data, it always has been a 
central question whether this agreement is merely the result of properly incorporating chiral 
symmetries, or if there is some deeper, underlying reason. One recent development that sug- 
gests that AdS/QCD is actually not that similar to real QCD is the conjecture of Strassler 
that event shapes in theories at large 't Hooft coupling are spherical [3]. This possibility was 
already strongly suggested by the result of Polchinski and Strassler that in hadrons at large 
't Hooft coupling, all partons are at small-x ("wee") [4]; thus a parton shower with highly 
energetic partons is very unlikely. The spherical-event conjecture was proven by Hofman and 
Maldacena, who in a pioneering work set out to understand the collider physics of conformal 
theories with a gravity dual [5]. What they found confirmed that these theories behave very 
differently from QCD: instead of jet-like events one finds energies that are distributed in a 
spherically symmetric way in the calorimeter. For related work see [6] and also [7]. One 
aim of this paper is to verify that the results of Hofman and Maldacena persist in the case 
of a theory with conformal symmetry broken in the IR as in the case of AdS/QCD. Indeed 
we find that introducing the IR cutoff does not substantially modify the basic results, and 
generic scattering processes in AdS / QCD-like theories will lead to events with large spheric- 
ities (as opposed to jet-like events in QCD). Of course this result is not the first qualitative 
disagreement between AdS/QCD and real QCD: the resonance spectrum of a simple extra 
dimensional model is also quite different from QCD, unless a soft-wall metric is assumed for 
the 5D background. We have checked that the spherical nature of the event shapes persists 
even when the simple RS background is replaced by the soft-wall background. These results 
suggest that there is a real disagreement between theories with a gravity dual and QCD, 
as one would have expected from the beginning since gravity duals are only calculable in 
theories with large 't Hooft coupling. 

In order to better understand what these differences can be let us be more specific and 
consider QCD as a function of the number of colors N and the 't Hooft coupling A = 
g 2 N. In particular, we are concerned mostly with the value of A at distances much shorter 
than the confinement scale, so that by "small A" we mean asymptotically free or at least 
asymptotically weak coupling (e.g. a theory with a Banks-Zaks fixed point in the UV). A 
brief survey of the (N, A) plane is illustrated in Fig. [I] At large N, and for any A, we expect 
a theory with some stringy reorganization in a 1/N expansion [8] (but not necessarily a 
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Figure 1: A survey of the N — X plane. 

gravity dual). On the other hand, for large A we expect the existence of a weakly coupled 
gravity dual, which is generically a string theory; however for large N the string modes 
decouple and a pure gravity dual should be enough to describe the theory. Thus we get the 
basic characterization of the four corners of the (N, A) plane: for large N, A the theory is 
reminiscent of a Randall-Sundrum [9] (or Klebanov-Strassler [10]) model: field theory (and 
classical gravity) on a warped background. For large A, small N one should find string theory 
on a warped background. For large N, small A we still expect a stringy description, but not 
necessarily higher dimensional, while for small N, A we get ordinary perturbative QCD, with 
its characteristic parton showers followed by hadronization. 

The main aim of this paper is to examine the various transitions in the (N, A) plane. At 
first sight changing N from small to large values looks like a very dramatic transition, since 
in the large N limit the decays of mesons will be 1/iV-suppressed. As a result, at strictly 
infinite N, the theory behaves vastly differently from a small N field theory. In infinite N 
QCD, an e + e~ collider would produce no hadrons at all, except when tuned to have precisely 
the energy of a resonance. However at small enough N, above the first few resonances, the 
spectral function will be dominated by a smooth, flat continuum. This is true at any A: the 
cross section at small N is continuum, but at large N is a discretuum. At very large but 
finite N, we can study production of narrow on-shell mesons which subsequently decay. In 
the large-A case these mesons are KK modes. In the small-A case, they are flux tubes. 

Arguments for the flux tube description of hadrons date back to the 1970s, and include 
in particular 't Hooft's characterization of phases of gauge theories based on the behavior of 
Wilson and 't Hooft operators [11]. Roughly speaking, there is a discrete set of possibilities 
for these behaviors, one of which is the confining area-law behavior of the Wilson loop, which 
suggests that the partition function has contributions of the form e~ A2LT . This suggests that 
there are states in the theory with constant energy density per unit length, so that their 
energy is E = A 2 L. Such states would look stringy. Various aspects of this argument have 



A — » oo 



A | 



A -» 



2 



2 



been made with various amounts of precision over many years; a partial list of references can 
be found in Ref. [12]. The upshot is that confinement, in the sense of area-law behavior over 
some range of distances, implies that hadrons are approximately described by the Nambu- 
Goto action. There are similar flux tube states at large 't Hooft coupling [13,14], but they 
are parametrically heavier, and we will not discuss them here. 

Our goal in this paper is to argue that the small N to large N transition, at given A, is 
fairly smooth and does not drastically alter qualitative properties of the physics, despite the 
singular nature of the infinite N limit; in other words, we expect to have a significant overlap 
between partonic, QCD-like descriptions, and hadronic, narrow-resonance descriptions. This 
supports the picture of "parton-hadron duality." On the other hand, the transition from 
small to large A is in some ways more mysterious. Large 't Hooft coupling produces deeply 
bound meson (and glueball) states that do not resemble the experimentally observed mesons 
(in particular, they do not fall on linear Regge trajectories), but which are calculable in 
supergravity duals [15-20]. The exotic properties of light, low-spin states at large 't Hooft 
coupling have been emphasized in Refs. [4,17,18]. Certain decay modes and form factors are 
strongly suppressed at large 't Hooft coupling [18,21]. Other aspects of the large 't Hooft 
coupling physics have been studied that contribute to the picture that strong coupling at 
short distances leads to novel physics [22-24]. The most recent such observation is that, 
whereas small A leads to the familiar jet structure in QCD-like events, large A leads to very 
spherical events as proved by Hofman and Maldacena. As far as we are aware, there is no 
tool currently available that can tell us whether the spectrum and event shapes transition 
smoothly as a function of A or change abruptly at some critical A. The latter possibility 
may seem strange, but a similar transition happens in the spectrum of modes of a static flux 
tube in Af = 4 SUSY Yang-Mills [25]. At large A many issues involving very high energy 
scattering of hadrons, such as Regge trajectories, the Pomeron, and saturation [4,26,27], 
have been studied. Despite substantial quantitative differences, and some qualitative ones 
(e.g. the Pomeron couples to individual partons at small A, but to hadrons at large A), many 
aspects of the physics of high energy scattering (e.g. the saturation of the Froissart bound) 
are similar at small and large A. Nonetheless, the details of dynamical questions at energies 
near the scale of the lightest resonances are strikingly different at small and large AQ We 
should emphasize that this idea occurs, in various forms, throughout the literature, and is 
not new to us. However, we hope to present particularly simple and clear illustrations of 
the different regimes and make it apparent that these differences are generic, not artifacts of 
particular examples that have been studied in the past. 

The paper is organized as follows: in section [2] we study the connections between the 
predictions of perturbation theory (PT) for a two-point function of currents, (J(x)J(O)), 
and the large N spectrum. At any A, the deep Euclidean behavior of such a two-point 
function is determined by conformal symmetry, and unitarity relates the Minkowski and 
Euclidean limits. Various possible spectra are consistent with these two facts. However, 
in the perturbative case, we expect the leading-order deep Euclidean calculation to capture 
the physics almost everywhere in the complex plane, with small corrections. We confirm 

1 Another argument against the existence of a gravity dual of QCD was recently made in [28]. 
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this expectation by showing that a large-iV Regge-type spectrum, rr? n ~ n, is consistent 
with PT in most of the complex plane, except very near the Minkowski axis. At finite N 
the poles become cuts, but the behavior of the two-point function in most of the complex 
plane will be similar. On the other hand, the RS-type spectrum, m n ~ n, deviates sharply 
from PT in much of the complex plane. There is nothing wrong with this, but it implies 
we should only expect such a spectrum in theories with large 't Hooft coupling, where PT 
breaks down very badly. The analysis of the spectrum continues in Section [3j which presents 
a simple quantum-mechanical toy model that illustrates how the spectrum of light states can 
change continuously from flux tube-type to RS-type as a coupling constant is varied. This 
suggests that the KK-like modes that are characteristic to AdS/QCD are actually not QCD 
states, but rather they are deeply bound mesons, that are more related to the short-distance 
Coulomb-like potential, and not to the long-distance confining linear potential. 

Having established the large- N spectrum looks quite different at large and small A, we 
turn our attention to more detailed properties of scattering events. We look at large but 
finite N, so that we can produce a narrow resonance on-shell and watch it decay in a sequence 
of steps. First we examine decays of Kaluza-Klein modes in an extra dimension in section 
|4| We show that the end result of the decay chain is a large number of the lightest modes 
moving in approximately independent directions, so that events appear spherical. We repeat 
the analysis for a soft-wall type background, and show that the spherical nature of the 
event shapes remains. Although we cannot reliably calculate the same process at small N, 
where stringy corrections would be important, we note that the physics is similar to that 
computed by Hofman and Maldacena, in which conformal theories at large N and A give 
rise to spherical events. We speculate that large sphericities are characteristic of any large 
A theory. 

In section |5j we repeat this analysis for a different model of hadrons, namely quasi-stable 
flux tubes as we would expect to find in small A, large N theories. We study models of flux 
tube evolution and breaking (motivated by Schwinger pair-production in the chromoelectric 
field), and find that the resulting events have a characteristic "jetty" structure, with energetic 
particles moving in opposite directions and only relatively soft particles in between. This 
is very similar to the well-known jet structure that arises from perturbative QCD and the 
parton shower, independently of any assumptions about hadronization. We take this as 
evidence that the small- N /large- N transition at small A is much smoother than the small 
to large A transition at large N. 

These pieces of evidence, taken together, provide support for our picture of the (N, A) 
plane. The old idea of parton-hadron duality is plausibly true everywhere in that plane, with 
the main effect of varying N being to broaden mesons until they merge with a continuum, 
with a large region of overlap between partonic and hadronic descriptions. The effects of 
raising A are much more dramatic, with the spectrum of light excitations and the shape 
of events in a detector being starkly different in the small and large A limits. We have 
no information about what the transition region is like; it may be smooth, with no phase 
transition, but it is at least a crossover between very different regimes. This suggests that 
if AdS/QCD is ever to be useful for truly QCD-like theories, we will need qualitatively new 
ideas and a better understanding of this transition region. 
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2 Perturbation Theory and the Large- N Spectrum 



We are interested in understanding qualitative properties of confining gauge theories in the 
plane of the number of colors N and the 't Hooft coupling A. Because we are interested 
in qualitative properties, we won't concern ourselves too much with which theory we are 
working in. In particular, for QCD it does not make sense to take A very large, unless 
one considers a brane construction of a QCD-like theory with new degrees of freedom as in 
Ref. [29]. For concreteness, one could imagine studying the M = 1* theory obtained by mass 
deformation of M = 4 super- Yang-Mills, which has confining vacua in which one can vary iV 
and A [30,31]. In any case, we never have a fully controlled calculation over the whole plane, 
so we will have to pick certain corners and study them. In this paper we will be content 
with simplified calculations that exhibit general qualitative properties of the corners of the 
(N, A) plane that we are interested in. 

In a theory with small 't Hooft coupling at short distances, we can do a reliable pertur- 
bative calculation of the two-point function of a current in the deep Euclidean region. We 
trust from the operator product expansion that corrections to the calculated behavior are 
controlled: they are either proportional to a small coupling, or they are power-suppressed. 
The total cross section of a process like e + e~ — > hadrons is computed by analytically con- 
tinuing the Euclidean answer to the Minkowski region and taking the discontinuity across a 
cut. At large N, we know that the cross section vanishes except exactly on a resonance. This 
tells us that the large- N behavior of the two-point function is controlled by a meromorphic 
function, i.e. one which is analytic throughout the complex plane except at a discrete (but 
infinite) set of poles. It turns out that there are different meromorphic functions that have 
the same deep Euclidean asymptotics. However, we will argue that some of these meromor- 
phic functions are much closer to the perturbative answer over the whole complex plane than 
the AdS/QCD-like meromorphic function, so already PT for the two-point function seems 
to be telling us something important about the spectrum of hadrons. 

Two decades before AdS/CFT or any examples of confinement with large 't Hooft cou- 
pling at short distances, Migdal attempted to find a large N ansatz for the QCD spectrum by 
Pade approximation of perturbative Euclidean calculations [32]. The idea behind the Pade 
approximation is to produce an analytic continuation of the perturbative result for H(Q 2 ) 
into the non-perturbative regime. The leading order perturbative result for the two-point 
function of two conserved currents is just given by (dropping a factor of for convenience) 



with s = Q 2 . In the Pade approximation one is looking for a ratio of two polynomials 
rip(s) = Pn(s)/Qn(s) which at a given point (for example s — 1) reproduces the first 
2N + 1 terms in the Taylor series of n PT . One can then take a limit as the order of these 
polynomials increases, finding a function with infinitely many poles. The resulting function 
obtained this way for iV — > oo is 




(2.1) 




My^/P) log(-s//3) - nY (^/P) 



(2.2) 
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Here f3 is chosen as {2x\)~ 1 where X\ is the smallest root of Jq(x); this is chosen such that 
Ilp(s) has its first pole at s = 1/2. We can see in Fig. |2]that it is a very good approximation 
to the PT result for large Euclidean momenta (by construction). 
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Figure 2: Perturbation theory vs. meromorphization at Euclidean momentum, U(Q 2 ): black 
solid line, logQ 2 ; green dashed line, Pade approximation; red dotted line, digamma function. 

It is well-known that this Pade approximation fails to capture some of the essential fea- 



tures of the non-perturbative physics of QCD. For example, the poles of Eq. (2.2) appear 
at vn? n ~ n 2 , in stark contrast with the expected Regge answer m 2 ~ n. The difference be- 
tween two-point functions with these types of behavior in the Minkowski region is illustrated 
in Figure [3j While Migdal's answer does not match the expected QCD-like behavior, it is 
remarkably precisely what one would calculate in the AdS/QCD framework [2], which can 
be considered as a simple model of what this quantity would look like in a theory at large 
't Hooft coupling. This correspondence was first pointed out by Shifman in Ref. [33], and 
has been further discussed in Refs. [34]. Thus, while Migdal failed to produce a good model 
of the QCD two-point function, he stumbled upon a good model of large 't Hooft coupling 
theories! 

As a contrasting model, one could consider the digamma function, ip(x), which is the 
logarithmic derivative of the gamma function, 

*M = £g = | (ftT " TTk) ' (2 - 3) 

This function has the asymptotic expansion: 

1 °° B 

^( :r + l)«tog(x) + _-£^jL (2.4) 

n=l 

so that it has the right sort of behavior to match perturbative asymptotics. This function 
has been discussed as a model of a QCD-like spectral function a number of times in the 
literature, beginning with Ref. [35]. It allows us to build a model of the spectral function 
with the first pole in the same location (setting the p mass), s = 1/2, as: 

U R (s) = V> (s - i J . (2.5) 
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Note that this function, unlike the Pade result, has poles evenly spaced in m 2 , not in m, so 
that it matches expectations for QCD. We show this function in the plots for comparison. 




Figure 3: Different mass spectra with similar deep Euclidean asymptotics for II(s): green 
dashed line, analytic continuation of the Pade approximation; red dotted line, digamma 
function. 




Figure 4: Smeared spectral function, ^(H(s + i) — — i)): black solid line, perturbation 
theory; green dashed line, Pade approximation; red dotted line, digamma function. The 
digamma function with proper Regge physics provides a reasonably good match, while the 
Pade approximation oscillates about the perturbative result with large amplitude. 

The obvious question is which of these approximations one should take more seriously 
for the case of large N and small 't Hooft coupling. On the face of it both of these functions 
reproduce the perturbation theory results quite precisely. However, we will argue that that 
despite the good agreement with PT for the deep Euclidean region the Pade approximation 
actually does not match PT over the whole complex plane all that well, and that other func- 
tions like Il#(s) exist that will capture the global features of PT much better. To understand 
why the Pade approximation fails we need to understand the region of momenta where the 
Pade approximation is bad. Clearly it is doing extremely well for the deep Euclidean region. 
However, some of the Minkowski region should also be attainable through PT, at least when 
we are away from the poles, i.e. away from the real axis s > 0. One particularly simple 



7 



quantity of this sort is a smeared version of the spectral function p: 

p A (s) = i (n(a + iA) - U(s - iA)) . (2.6) 

A function of this sort has been discussed by Poggio, Quinn, and Weinberg [36]. Many 
comparisons of experimental data on QCD with PT implicitly rely on an approximate local 
quark-hadron duality that allows a result smeared over hadronic resonances to match PT. 
At low orders of PT, the choice A ~ 1 is effectively an infrared cutoff at about the scale 
Aqcd- (In fact, since s = 1/2 here corresponds to the p mass squared, this acts like an 
IR cutoff of about 1 GeV.) However, the Pade approximant disagrees quite sharply with 
the expected PT answer, as can be seen in fig. |4} It would be rather surprising if such a 
smeared perturbative calculation could disagree so strongly with the true hadronic answer 
at large values of Minkowski momentum! On the other hand, the smeared spectral function 
calculated from the digamma function Ur will still be a very good approximation to the 
PT answer even for Minkowski momenta. The agreement is shown over the complex plane 
in Figure [5j Note that 11^ has larger error in a A^cn-size region around the origin, but 
otherwise disagrees only in a narrow strip along the Minkowski axis. (Of course, one should 
not expect that the agreement is good until \Q 2 \ ^> Aq C d, so this is entirely as expected!) 
On the other hand, the Pade result lip has large error in a wedge that grows progressively 
larger as one goes deeper into the Minkowski region. 

One would hope that one could complete Migdal's program of uniquely fixing the appro- 
priate meromorphic function relevant for large N QCD by also requiring a good match to 
the Minkowski region. An attempt in this direction is described in detail in Appendix A, 
where we do show that the spacing between the poles rn\ +1 — m\ should not be growing, and 
that the simplest meromorphic function satisfying all constraints from PT has asymptotic 
properties agreeing with those of the digamma function. However, we could not exclude 
some more exotic possibilities, so we conclude that Migdal's hope of finding a systematic 
procedure to convert perturbative results to meromorphic functions does not seem tractable, 
at least when only the two-point function is considered. 

To summarize this section: we found that in order to apply Migdal's program of analyt- 
ically continuing the perturbative result one should use the digamma function rather than 
the Pade approximation. The use of such meromorphic functions with QCD-like spectra has 
been previously discussed; see, e.g., Refs. [35,37]. The result obtained by Pade approxima- 
tion is a good guide to the large 't Hooft coupling spectrum (that is the case when AdS is 
a good description). Of course, at large 't Hooft coupling PT breaks down, so it should not 
be surprising to see that there can be a qualitatively different spectrum for a large 't Hooft 
coupling. In the next section we will discuss a very simple toy model that illustrates how 
the spectrum can change as the 't Hooft coupling changes. 
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Figure 5: Contours of constant fractional difference, iog"( g °) > between the leading-order 
perturbative result and the full two-point function. At left is the Pade case lip and at right 
the digamma function lip. White regions correspond to large fractional disagreements, and 
dark regions to very good agreement. The digamma function disagrees with PT in a circle 
of order Aqqd, and very close to the Minkowski axis, while the Pade approximant disagrees 
in a growing wedge in the Minkowski region. 

3 A toy model of bound states: where have all the KK 
modes gone? 

We have seen that a QCD-like spectrum is much more consistent with perturbative calcu- 
lations than a Pade or AdS/QCD-like spectrum. It would be interesting to understand in 
more detail how the spectrum can change as one varies the 't Hooft coupling. In this section 
we will discuss a very simple toy quantum-mechanical model of bound states in which the 
spectrum can change between QCD-like and AdS/QCD-like. We make no claim that this is 
an accurate model of the dynamics of any quantum field theory, but it does provide some 
intuition about how strong coupling at short distances can change the spectrum. 

For studying the spectrum of mesons in a QCD-like theory, we expect that we have a 
confinement potential well-approximated by the Cornell form, 

V(r) = -- + ar, (3.1) 

r 

with cr the string tension. The spectrum of mesons is, at least qualitatively, determined 
by a relativistic bound-state calculation involving this sort of potential. (There have been 
various studies of meson and glueball spectra in Coulomb gauge QCD, which resemble these 
potential models, including some which are quite sophisticated [38].) 

What is notable about the large A examples is that the non-QCD-like states have masses 
of order M where M ~ X^^y/a <C \fo. This suggests that these bound states involve 
physics at distances where the ar term of the potential is not yet dominant. Rather these 
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states seem to be related to the short-distance part of the potential, and so they are likely 
better characterized as deeply bound mesons due to the Coulomb nature of the potential at 
short distances. In QCD, there is only one scale, Aqqd, and it's difficult to see how such 
states could arise. So let's turn our attention away from QCD and instead think about 
Af = 4 super- Yang-Mills deformed by mass terms to an Af = 1 confining theory. The mass 
terms give one scale, m, at which we expect strongly-bound mesons at large A, where the 
string tension is \f\rn 2 . 

The static potential is computed along the lines discussed by [39]. In particular for RS 
(as discussed in [41]) it will be Coulomb-like until the string is long enough to reach the IR 
brane; at that point, any additional string length lies along the IR brane and the potential 
grows linearly. In other words, 

V(r) = { ~f v r ^ rc (3.2) 
v ' \-f +cr(r -r c ), r > r c v ' 



Here 

Air 2 (2g 2 YM N) 1/2 



c = 



rQ) 4 



0.3231^, (3.3) 



as in the AdS / CFT calculation. We can read off from Maldacena's calculation of the static 
quark- ant iquark potential that the string just touches the IR brane at z IR when the quark 
and antiquark are separated by a distance 

r c = ZlR ^J\2 ~ (3.4) 

Finally, the string tension is given in terms of the 't Hooft coupling and the location of the 
IR brane: 

i R 2 ( 

° 2na'z 2 IR 2nz 2 IR 1 ' 

Our approach will be to use this potential to model the confining force between a scalar 
quark and anti-quark in a simple toy model of relativistic quantum mechanics. That is, we 
solve for radial excitations that are eigenstates of the Hamiltonian 

H=y/p* + m? + V(r), (3.6) 

where p r is the radial component of the momentum. One can think of this as a very rough 
approximation to a Hamiltonian formalism for the theory in Coulomb gauge, approximating 
the Coulomb potential by the Wilson-loop potential, although strictly speaking they are 
known to differ [40]. 

We will make one further approximation in order to do more of the problem analytically, 
which facilitates an easy numerical solution for the eigenstates. Namely, we will use a 
variational ansatz with basis functions of the form e _/3r L 2 (2/3r), and we would like a potential 
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Figure 6: Comparison of the RS potential (red, dashed), shifted by a constant value of 10.5, 
and our analytically useful replacement of it (blue, solid). 

V(r) that has analytically calculable overlap with these basis functions. In particular we 
take the parameters of the model to be specified by m — 1, a — 10, c = 3 (note that with 
these parameters y/X ~ 9.28 is large), and we approximate V(r) by the more analytically 
tractable: 

V mod {r) = — exp(-3r) + lOr, (3.7) 
r 

in effect changing the sudden ^-function like cutoff to a smoother exponential shut-off of the 
Coulomb part of the potential. 

This modified potential agrees quite well with the original RS potential, shifted by a 
constant value of 10.5, as illustrated in Figure |6j Because there is no clear correct definition 
of the zero of the potential, we will ignore the constant shift. (Of course, if we look for evenly 
spaced masses, a constant shift in energy eigenvalues is irrelevant; if we look for evenly spaced 
squared masses, a constant shift will matter, but will become less important at high enough 
excitation number.) 

Finally, to get a toy model of the transition from a QCD-like spectrum to an RS-like one, 
we vary the coefficient of the first term in V mo d{r): 

V 7 (r) = -- exp(-3r) + lOr. (3.8) 

At 7 = we have only a linear term, resembling a QCD-like theory where Coulomb at- 
traction plays very little role in the pattern of bound states, while at 7 = 3 we have a 
theory which more closely resembles RS. The squares of the eigenvalues of the Hamiltonian 
(r) are plotted for 7 = 0, 1, 2, 3 in Figure [Tj One can see that as the strong 
Coulomb attraction is turned on, the low eigenvalues become closely spaced and take on a 
qualitatively new character (although they change continuously). These are deeply bound 
states, sensitive to the strong short- distance Coulomb attraction and relatively insensitive 
to the long-distance linear confinement. They are the toy model analogs of KK modes. 
Higher eigenvalues are spaced successively farther apart, eventually merging into the origi- 
nal pattern of linearly spaced mass-squared. They are the toy model analogs of resonances 
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on linear Regge trajectories. Note that at 7 = one finds roughly even spacing consistent 
with m 2 (n + 1) — m 2 {n) = — . 
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Figure 7: Spectrum of our toy Hamiltonian . 

A more reliable and field-theoretic approach to modeling the effect of strong short- 
distance Coulomb attraction on bound states would be to use Bethe-Salpeter equations, 
but already our toy model has interesting physics that suggest a picture of how the transi- 
tion from QCD-like to RS-like dynamics can happen. As a final remark on this topic, we 
should point out that literature beginning from a "soft wall metric" background [42] com- 
putes AdS/QCD potentials which are good matches to lattice data and which work well 
phenomenologically in Salpeter calculations [43] . On the other hand, we will see in the next 
section that event shapes in soft wall backgrounds are not QCD-like, so if such potentials 
are to fit into a coherent framework for QCD phenomenology, one would have to understand 
stringy corrections to the effective action. 



4 Dynamics for large N, large A: Spherical Event Shapes 

So far we have been discussing the two-point function of a current with itself, which de- 
termines the spectrum of associated resonances and the total cross section of a process like 
e + e~ — > hadrons. However, the real physics of e + e~ — > hadrons is much richer: in QCD, for 
instance, we know that such an event will lead to a parton shower that produces collimated 
jets of particles moving in a small number of directions. We now want to turn our attention 
away from the crude questions of spectrum and cross section and toward the richer ques- 
tions about the structure of events. This will mean going to finite N and considering how 
resonances decay and broaden into continuum structure. 

Our tool for understanding whether an event is spherical or jetty will be the "sphericity" 
[44], which is defined as S = §(A 2 + A3) where A 2 and A 3 are the two smallest eigenvalues of 
the matrix: 
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where p( 1 ' is the 3-momentum of the 2 th final-state particle. The sphericity is 1 for a com- 
pletely spherical event, and for two back-to-back particles. Perturbative QCD predicts 
small sphericities of order fewxO.Ol. 

Extra dimensions, in general, predict a sizable sphericity. To see why, recall that in a 
flat periodic extra dimension, or a flat interval with Neumann boundary conditions on both 
ends, there is a discrete conserved 5D momentum, or "KK number." The wave functions 
are cosines and triple overlaps vanish except when the mass of the highest state is precisely 
the sum of the masses of the two lower states. This case is not very interesting; there is 
no phase space for decays. In a UED model, the exact KK number conservation is broken 
by boundary masses or by localized boundary kinetic terms. However, in this case, as in 
many other theories with an extra dimension, there is still an approximately conserved KK 
number. In particular, in a two-body decay of a heavy KK mode, the masses of the two 
daughters are very close to being equal to the mass of the parent; as a result, there is little 
phase-space and the daughters are not very boosted. Thus, if they subsequently decay, the 
new daughters will again not have much of a boost, and the momenta of the various particles 
in the final state will not be strongly correlated with any particular direction. This leads to 
spherical events. 

In particular, Figure [8] shows one example of the shape of such an event. This is the decay 
of a heavy KK mode of a gauge boson in flat space with a Neumann boundary condition on 
one end of the space and a Dirichlet boundary condition on the other. The spikes radiated 
out represent the momenta of stable daughters at the end of the decay chain, with length 
proportional to energy. Note that daughters are flying out in all directions, and if we were 
to place a spherical detector of large radius around the origin, it would see a fairly uniform 
distribution of energy. In particular, the sphericity is high. This is not a jetty event, but one 
closer to the uniform energy distribution found in the conformal case [5]. The distribution 
of sphericity for such decays is also shown in Figure [HJ where we see that an overwhelming 
fraction of events have fairly large sphericities > 0.1. 

Although approximate conservation of KK number is most obvious in flat space, it turns 
out that even in warped space it is still true. As a first example we consider a scalar field in 
an RS background with cubic interaction: 

S = j d'x j Z " dz W^j ( v ^d^d u 4>-d z 4>d z cf>)+(^j ss0 3 j, (4.2) 

with R the AdS curvature radius. We take Zjjv ~^ 0, requiring that wave functions be nor- 
malizable at the boundary, and we impose a finite zir where Dirichlet boundary conditions 
are imposed. This means that our wave functions are: 

Mz) = Z2 ' HmnZ) (4-3) 

^/-R 3 z] R J 1 (m n z IR )J 3 (m n z IR )/2 

where J2,{rn n ZiR) = and the denominator serves to normalize the wave function. 

From these wave functions, we can compute couplings as triple overlaps and then compute 
decay chains. While it is not obvious from inspecting the wave function, it turns out that 
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Figure 8: On the left: typical spherical event from the decay of the 200th gauge KK mode 
in flat space. On the right: the sphericity distribution from 6000 decays of the 200th gauge 
KK mode in flat space. 



the triple overlaps still show a strong tendency to almost conserve KK number, so that 
decays with little phase space are preferred. As an illustration of this, we plot in Figure [9] 
the distribution of the speeds, t>i ;2 = |pi,2|/-£a,2, of the two KK modes originating from the 
decay of the 100th scalar KK mode in RS. Notice that in almost every case, the decay is 
to two relatively slow-moving, unboosted particles; the masses of the two daughters sum to 
nearly the mass of the parent. Just as in flat space, the tendency to prefer decays with little 




Figure 9: Distribution of the speeds v\ and of the two Kaluza-Klein modes originating 
from the decay of the 100th scalar KK mode in RS. Note that typically the decay is to 
slow-moving particles. 



phase space will lead to very spherical events. A typical such event for the decay of a heavy 



gauge KK mode in RS is shown on the left panel in Figure 10, and the sphericity distribution 
for many such events is shown in the right panel in Figure 10 Because the zeroes of Bessel 



functions are very nearly linearly spaced, it turns out that the number of the KK mode is 
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almost conserved. This is illustrated in Figure 11 The 200th KK mode decays most often 
to two modes with mode number summing to 198, then to 199, then to 197, and so on. 




s 



Figure 10: Left: a typical spherical event from the decay of the 200th KK gauge KK mode 
in AdS. Right: Sphericity distribution event from 9000 decays of the 200th gauge KK mode 
in AdS space. 



Figure 11: Branching fractions for the 200th KK mode with the KK number of one of the 
daughters on the horizontal axis. In the top line (Blue) the KK numbers of the daughters 
add up to 198, in the second line (Red) they add to 199, in the third line (Green) they add 
to 197. 

RS backgrounds are similar to flat space, to the extent that they predict mass spectra 
with nearly linear spacing. What happens if we tune an extra-dimensional theory to achieve 
the Regge spacing w? n ~ n, which holds in small 't Hooft coupling theories as discussed in the 
previous sections? One might hope that almost-conservation of KK number is lost and event 
shapes are very different. We will now see that is not the case; despite the very different 
spectrum, decays with little phase space are still preferred and sphericities are large. Our 
starting point will be a gauge field in the by now familiar "soft-wall" model of AdS/QCD, 
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in which the metric is AdS but the action is [45]: 



S 



d 5 x e -^-^F MN F MN , 



-A 2 z 2 



(4.4) 
(4.5) 



The exponentially damped dilaton leads to an equation of motion that can be rewritten as a 



Schrodinger equation for the eigenvalues vnr n with harmonic oscillator potential 



~ z' 



at large 



z, which yields the desired spectrum. In the original soft-wall model, the dilaton is simply 
assumed to have this profile, due to some unknown dynamics. 

In this model, triple overlaps of wave functions vanish; we must deform it slightly to get a 
model in which we can calculate decays. As a simple deformation, we include a "UV brane" 
at z = zuv and impose Dirichlet boundary conditions there (the other boundary condition is 
simply normalizability at z — > oo). This gives a small deformation to the wave functions and 
masses and allows us to compute decay chains. It turns out that the correct wave functions 
are given by: 



^n{z) = C n Z 2 U 



~™>r. 

4A 2 



+ 1,2,*' 



j z uv 



(4.6) 



de- 



with U a confluent hypergeometric function. The condition U yj^r + 1,2 

termines the masses m n and the constant c n is chosen to give canonically normalized KK 
modes. 




Figure 12: At left: a typical spherical event from the decay of a heavy (the 46th mode) 
soft-wall gauge boson. The UV brane was chosen to be at zuv = 0.01/Aqct>. At right: 
average sphericity in 500 decays of the n th KK gauge boson in the soft-wall background. The 
horizontal axis is the KK number and the vertical the average spehricity of 500 decays of 
that mode. 



In Figure 12 we show a typical spherical event, and the distribution of sphericities. The 
numerical treatment of the hypergeometric [/-function is much harder than that of Bessel 
functions, and due to the smaller spacing among the masses there are fewer kinematically 
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allowed decays for a given KK mode. Nevertheless, even with a relatively low number of 
events the trend toward final states with larger sphericities is quite clear. This is again 
because despite the very different mass spectrum, the wavefunction overlaps prefer the de- 
cays with small phase space and correspondingly small boosts. The soft wall predicts more 
spherical events which are un-QCD-like, despite its QCD-like mass spectrum. Attempts to 
put the soft-wall background on a dynamical footing, as in Ref . [46] , seem unlikely to remedy 
this problem unless they involve inherently stringy dynamics. One perhaps more promising 
direction, begun in Ref. [47], is to begin with the solvable large- N 't Hooft model in 1+1 di- 
mensions and attempt to reformulate it as a 2+1 dimensional theory. Such an approach could 
also clarify the relationship between holographic wavefunctions and the internal structure of 
hadrons [48] . It remains to be seen how far this program can be pushed, and whether it will 
lead to anything resembling the soft wall. Also, note that demanding "asymptotic freedom" 
in the sense of a 5D dilaton which approaches a weak coupling value will not change these 
event shapes. To truly model the dual of a weak coupling theory, one would have to include 
light fields (string modes) dual to the many operators with small anomalous dimension. For 
now, it seems that there is still a major obstacle to obtaining QCD-like dynamics from a 
5D theory, one which will likely require new techniques in string theory to be overcome. 
For now, we will turn back to 4 dimensions, and look at how much older models of hadrons 
generate event shapes that are much closer to those of QCD. 

5 Large AT, Small 't Hooft Coupling: Flux tubes 
5.1 Flux tube Breaking, Directionality, and Jets 

If, as we have argued, the mesons of RS do not smoothly match onto the mesons of a 
QCD-like theory, what sort of mesons do? There is a long history of modeling the highly 
excited states of QCD-like theories using flux tubes. Among the reasons for doing this are 
the idea that large N QCD is a string theory because it has a topological expansion [8], the 
experimental fact of Regge trajectories [49], the emergence of flux tubes on the lattice [50], 
and general arguments about the behavior of 't Hooft versus Wilson loops [11, 12]. Unlike in 
theories with large 't Hooft coupling and good gravity duals, the stringy physics in confining 
theories at small 't Hooft coupling begins at the scale of the lightest resonances; there is no 
large y/X to provide a parametric separation of scales. 

Before we look at more detailed modeling of flux tubes, let's begin by thinking about the 
first breaking of the flux tube that would be produced, for example, in an e + e _ collider. We 
begin with a quark and antiquark moving in opposite directions, each with momentum p, 
and assume they feel a linear potential, i.e. a constant attractive force a toward each other. 
We assume the probability of breaking is uniform along the string length; thus, as the string 
grows linearly, initially the decay probability per unit time goes as t 2 . If the string breaks 
immediately, then there are massless quarks moving in opposite directions, and we have a 
clean two-jet event, just as in PT. If the string breaks when fully extended, we have two 
heavy static flux tubes of mass 2px and 2p(l — x), with x the fractional distance along the 
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Figure 13: Distribution of the speeds v\ and Vi of the two flux tubes originating from 
breaking of an initial flux tube. The state begins with back-to-back quark and antiquark of 
15 GeV momentum, and the string tension is 0.177 GeV 2 . The decay rate per unit length 
is, from left to right, 0.01 GeV 2 , 0.0025 GeV 2 , and 0.000625 GeV 2 . This corresponds to 99% 
of the strings having decayed by the time they are one-quarter, half, and fully extended, 
respectively. Note that even in the latter case many decays are to fast particles. 



string where the break occurs. These flux tubes will subsequently decay, but the event will 
not be jetty. After the time r osc — p/a at which the string is fully extended, the quark and 
antiquark turn around and again pick up speed, so at longer decay times there can still be 
many events at which the decay products are boosted. 



In Figure [13] we show the distribution of speeds of the decay products as the decay rate 
decreases. Provided the string decays in a time smaller than that needed for the string to 
fully extend, the first breaking is to high-speed objects, and we expect that their further 
decays will widen the jets but maintain the qualitative "jettiness" of the event. The thrust 
axis of the jets will be the direction along which the flux tube was stretched, established by 
the initial momentum of the quark and antiquark, just as in PT. Note that in order for the 
string to decay before the quark and antiquark reverse direction, we need Td ec < t osc . The 
decay time will be r^ ec ~ , so this tells us that we need A" less than about This is 
consistent with large A" provided we look at strings that are energetic enough. 

If the decay rate is sufficiently small, so that T& ec > t osc , then we have slow daughters 
much more often, although the distribution is still very different than that for KK modes 
in Figure [9] Note, however, that for such small decay rates, the discreteness of the large A" 
spectrum invalidates our way of thinking about the decay. Given m 2 ~ nA 2 , we have 

A a . . 

m n+1 -m n w — -= ~ — . (5.1) 

>n m„ 



Thus the spacing between subsequent modes is of order the inverse of the oscillation time 
of our semiclassical string; if r^ ec > t osc , the width is smaller than the interval between 
subsequent modes, and it no longer makes sense to assume a continuous flux tube that can 
break at any location. In this situation a quantum mechanical treatment is required, and 
the resulting discreteness of the spectrum will play a role in calculating the decays. We are 
unable to treat such long decay times in the context of this semi-classical model. 
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5.2 A Simple Model of Flux Tubes 

We have simulated a very simple toy model of the decay of an excited flux tube, in order to 
plot some event-shape distributions. The literature has many more elaborate models, with 
varying degrees of input from PT. Some early references are [51,53-55]; more recently, the 
"gluon chain model" [56] has provided an example of a model that reflects both partonic and 
flux tube limits of the theory, and is consistent with important nonperturbative features of 
Yang-Mills theory like Casimir scaling at intermediate distances and iV-ality scaling at long 
distances. (For a recent review of confinement in QCD-like theories, with an emphasis on 
the role of these features, see Ref. [57].) 

Our model is a semiclassical one, based on linear confinement and decay by Schwinger 
pair production, following Ref. [51] as amended by Ref. [52] so that energy is conserved when 
the flux tube breaks. We begin with a quark and antiquark moving at high momentum in 
opposite directions along the z axis, separated by a constant chromoelectric field filling a 
cylinder of radius R, so that they are attracted by a constant force a (the string tension). 
The total energy of the system is (taking the quarks to be massless): 

E = ^2\Pq\+^2\Pq\+^2 <r Wq-Xq\, (5.2) 

where the last term is a sum of the energy from string tension over quark-ant iquark pairs. We 
do not reconnect pairs, so that each quark is associated with precisely one antiquark partner 
at all times (changing only when there is a decay). The flux tube can decay by the production 
of a quark-antiquark pair at some point along the axis between quark and antiquark. The 
newly created pair will have transverse momentum px (orthogonal to the flux tube axis) and 
the quark and antiquark will be separated by a distance d(px) = 2Et/&. The stretch of string 
that has been annihilated between the two quarks then compensates for the energy added 
by their initial momentum. Transverse momenta are distributed according to a Gaussian, 

cLN(pt) ~ d 2 px exp ^— -^V Our simulation takes small time steps, alternating a classical 
evolution step and a random decay step. In the classical evolution step, quark-antiquark 
pairs evolve according to their current momenta and the color force between them, treating 



Eq. (5.2) as a Hamiltonian governing the evolution. (We never reconnect pairs, so each 
quark remains connected to the same antiquark until the string fragment splits, generating 
two shorter strings.) The random decay step decays each string fragment with probability 
proportional to the time step and to the length of the string fragment; decays happen 
with px (orthogonal to the quark-antiquark axis) chosen from the Gaussian distribution 
and location of the decay chosen from a uniform distribution on the string consistent with 
the newly produced pair having the correct separation d(pr) and remaining between the 
original quark and antiquark. This is a highly simplified model, neglecting for instance any 
internal dynamics of the string, but we expect it to get the grossest qualitative features 
right, at a level that allows an approximate comparison with the RS results and with QCD. 
It conserves energy (because the newly created pair are always separated by an appropriate 
distance along the string axis) and momentum (because the newly created pair have equal 
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and opposite momentum), but is not fully relativistically invariant^ In Section 5.3 we will 
discuss the relativistic string action and show that the first few decays computed from it are 
qualitatively consistent with the results of this simpler model. 

For convenience we use the same string tension as in Ref. [51], o = 0.177 GeV 2 . We 
simulated a flux tube initially consisting of a quark and antiquark produced at the origin 
moving with p z = ±15GeV along the z-axis. The overall decay probability is suppressed 
from the original calculations to take into account the effect of large N: for the plots shown 
the probability per unit time and unit length of a break is = O.OlGeV 2 . We take time 
steps of 0.025 GeV -1 . As a stopping criterion, after each decay we flag as stable any new 
string segment that has 



m 



y/E(\xi - x i+1 \,Pi,p i+1 ) 2 - | pi + Pi+i| 2 < 1.54 GeV. 



(5.3) 



As we can see in Fig. [14] these events are quite jetty, with very low characteristic sphericities, 
as expected in QCD. 




Figure 14: Left: a typical flux tube event from the decay of a 30 GeV flux tube in the toy 
model described in the text. Right: Sphericity distribution from 5000 decays of a 30 GeV 
flux tube in the toy model described in the text. 



5.3 Including String Dynamics 

We can move beyond this crude model of flux tubes by including stringy dynamics: we will 
model the flux tube worldsheet with the classical Nambu-Goto action. We specialize to static 
gauge, t = r, and parametrize the spatial coordinate s with the conditions 

d s X-d T X = 0, (5.4) 
(d s xy+ (d T Xy = 1, (5.5) 
which reduce the string action to 

S NG = ^f drds (l + (d s Xy - (d r X) J , (5.6) 

-'In Appendix B we show that for large enough center of mass energies it does give a Lorent invariant 
description. 
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so that the equations of motion are simple wave equations (for a detailed discussion, see the 
textbook [58]). Note that with this parametrization, the string energy is simply the tension 
a times the length of the s interval. 

There is a subtlety here: if we use just the Nambu-Goto action for the string, with 
Neumann boundary conditions, the string endpoints will always move orthogonal to the 
string direction. However, we begin with a string stretched along the z-axis between two 
quarks with momentum in the z-direction. The Nambu-Goto action alone will not suffice 
to describe this system. This puzzle is resolved by noting that the quarks on the ends of 
the string are physical particles, and we should add the action of a massive point particle 
at each string endpoint to the Nambu-Goto action, with the constraint that the string must 
end at the location of the point particle]^] There is a recent discussion of this in Ref. [59] in 
the context of "quirks" , hypothetical long strings of some new gauge group with Mq 3> A, 
which we will follow for the key points. As in our simple Hamiltonian toy model, the quark 
(and antiquark) will experience a force, but now we can have a proper relativistic and local 
description: the force is related to the direction the quark moves in, and the direction tangent 
to the string at its endpoint. Specifically: 



-HaA-^ + ^L^I, (5.7) 



where v» is the component of quark velocity along the string direction and v± is the orthoj 
onal component. The corresponding rate of change of the quark energy is 

^ = -a (5.8 
dt x/l^ 



If this force is acting to slow down the quark, then energy conservation will require that 
the quark leaves a new bit of flux tube in its wake. If the quark is speeding up, then it is 
eating some existing flux tube as it moves. Thus adding the point particle action at the 
string boundary is, from a certain point of view, an extreme deformation of the Nambu- 
Goto dynamics; we're not just studying a classical string that is moving relativistically, 
we're studying a string that grows and shrinks as it moves relativistically. As noted above, 
the static gauge parametrization that makes the string equation of motion simplest has the 
feature that the length of the s interval is the string energy in string units. To continue 
to use that parametrization in this context, we have to allow the size of the s interval to 
grow or shrink: if the string energy, in units of the tension, is e(t) = E str ing(t)/a, then we 
can use a parametrization in which, for the initial back-to-back string, s takes values in 



2 Note that in the AdS/CFT context, one deals with strings that are described by a supersymmetric string 
action with D-brane boundary conditions moving in some curved space. In particular, there is no added 
point-particle action of the type we are considering here. However, AdS/CFT can still describe open strings 
that are dual to a quark and antiquark moving back-to-back with a flux tube forming in between. The 
resolution of the puzzle in this context is that the AdS/CFT dual involves a string that is falling in the 
AdS radial (holographic) direction. For a more detailed discussion, including some remarks about how AdS 
radial evolution can be dual to the parton shower, see Ref. [3]. 
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[— e(t)/2, e(t)/2]. More generally, after a string breaking, the quark and antiquark will lose 
energy at different rates, so s ranges in [— £i(t), £2(t)], with 



dei 

~dt 
de 2 
~dt 



1 dE„ 



a dt 
1 dE, 



a dt 



(5.9) 
(5.10) 



We will denote £i(t) + £2^) simply by e(t), which measures the string energy in units of o. 

Numerically, it's more convenient to use a coordinate living in a fixed interval, s = 
^^(t)^ ~ 1 ^ ( — ^' Because this new coordinate is a function of both the original s and 
t = t, the action will become more complicated: what was the derivative with respect to r 
at fixed s, for instance, is now a combination of derivatives with respect to s and r: 



dX 



dX 
ds 



ds dX 
dt dr 



ds 
dt 



des + 1 
' dt e 



2de l 
e dt 



(5.11) 
(5.12) 



Note that % is interpolating between the rate of quark and antiquark energy change, 
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dX 
ds 



ds 
dt 



accounts 



appropriate units, as we move along the string. Also note that the term 

for the fact that the quark can have a velocity component vu along the string direction, even 
though the string itself can carry only transverse momentum. 

Thus we can rewrite the action (5.6) in the (s,t) 
as it involves E\ 2 (t) 



coordinate frame. It's somewhat ugly, 



and their time-derivatives, so in this parametrization there appears to 
be some nonlocal dependence on what is happening at the string endpoints. Nonetheless, the 
key point is that the change in string energy is entirely determined by the change in energy of 



the quark and antiquark, as in Eq. (5.8). Thus, numerically, we have a simple procedure to 



follow. We keep track of a fixed number N^t of "string bits," i.e. locations at fixed intervals 



in s along the string. In a small time interval dt, we can first calculate from Eqns. (5.7) and 



(5.8 ) the change in momentum and energy of the quark and antiquark at the ends of a string. 



This then tells us how e{t) is changing. We then use the equations of motion computed from 
the action Eq. (5.6), written in the (s, t) frame and discretized, to evolve each "string bit" 
forward dt in time. In practice, we do this time evolution with a fourth-order Runge-Kutta 
algorithm. We then make small adjustments to the positions X(s) and velocities d T X(s) of 
the string bits so that the parametrization conditions, Eqns. (5.4) and (5.5), are satisfied 



as precisely as possible numerically, because in practice small discretization errors in these 
conditions seem more numerically troublesome than discretization errors in the equations 
of motion themselves. This gives us a slightly crude, but fully relativistic, algorithm for 
evolving a configuration of strings forward in time. 

The next step is to decay the strings. Here we use an ad hoc modification of the pro- 
cedure from Section 5J2 We continue to be guided by the Schwinger calculation, but 
we are now taking it somewhat further from its domain of validity. Namely, we assume 
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that at a given instant in time, the decay probability of the string is proportional to the 
length of the string, and that the decay is equally likely to happen anywhere in the coordi- 
nate s. Rather than transverse momenta, we pull an energy from a Gaussian distribution, 

dN(E) ~ d 2 E exp ^— ^-J, insisting that E > 2rn q , and remove a length of s interval that 

balances this energy. In general this interval may also have some momentum. To construct 
the momenta of the new quark and antiquark, we boost to the rest frame of the removed 
string interval and construct the momenta of two daughters of mass m q , perpendicular to 
the string direction at the center of the removed interval. (This ensures that in the case of 
the first breaking, this algorithm matches the simple semiclassical flux tube algorithm of our 
previous discussion.) Boosting back, we compute the components vu and v± of the velocities 
of the new q and q and the corresponding rate of energy loss. We divide each daughter string 
into Nut bits, just as the parent string had, interpolating between bits of the parent string to 
obtain the positions and velocities of the new bits. We place the final quark and antiquark 



at a location such that the constraint (5.5) is satisfied. 



Figure 15: Sphericity distribution from the string model, after the first four string decays 



The numerics of this procedure prove to be tricky, and we present results only for the 
first four decay steps. The results of this model are similar to those of the simpler static 
flux tube model. We have fixed Nut — 101 for the number of points we track along the 
string, and have taken time steps of size 0.003125 GeV -1 . The quark mass is taken to be 
0.3 GeV. Other parameters (tension, initial p z , decay probability per unit time and unit 



length in string units) are as in the previous subsection. In Figure 15 we plot the sphericity 
distribution. It is peaked toward low values, showing that the string simulation predicts 
highly directional, jetty events just as our simpler flux tube model did. 



6 Conclusions 

We have examined various aspects of large N gauge theories for small and large 't Hooft 
coupling A. We found that for small 't Hooft coupling agreement with perturbation theory 
suggests that the mass spectrum should asymptotically show the characteristic Regge-type 
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behavior rather than the KK-type scaling, that arises in the large A limit. We have given 
a simple toy model that sheds light on how one can interpolate between these two spectra. 
This model suggests that the KK-like modes are more likely to be deeply bound mesons 
in the Coulombic short-distance part of the potential rather than genuine QCD-like modes 
related to the confining long distance potential. The event shapes resulting from scatterings 
in the small and large A theories also seem to be qualitatively different. While for small A 
one expects QCD-like events forming jets, for large A we have shown that the final states are 
much more spherical, in agreement with the conjecture of Strassler and the calculations of 
Hofman and Maldacena for events in gravity duals of CFT's. We conclude that in order to 
fully capture the dynamics of QCD with small A inclusion of some sort of stringy dynamics 
seems to be unavoidable. 
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Appendix 

A Towards the completion of Migdal's program 

In this appendix we describe an attempt at completing Migdal's program of analytically 
continuing the PT results for QCD into meromorphic functions valid for all energy regimes, 
along the lines of discussion in Sec. [2} This will lead us toward functions whose pole structure 
seems to asymptotically agree with that of the digamma function, however as we will see 
some other more exotic possibilities can not be completely ruled out based on our arguments. 

Let us try to find a function that is a good fit to PT both in the deep Euclidean regime 
and the Minkowski regime away from the real axis (eg. also matches p& for A ~ 1). We will 
be looking for a meromorphic function of the form 



We require that this function (just as the Pade approximant) reproduces logs for — s ^> 1 
with no 1/s corrections. But in addition we also require that it be a good fit to pa(s) for 
s > 1. Let us first understand why the Pade approximation fails this requirement. The 
problem with the Pade approximation is that the succession of poles and residues obey 



oo 




(A.l) 
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rn: 



n 2 and also r n ~ n. The expression for pa is 



' (s — m^J z + A 2 

Clearly, for s = one gets a contribution of r n , while the other terms are highly suppressed 
by the denominator. If \r n \ is not asymptotically approaching a constant (or decreasing), 
then there is no way to correctly approximate for large n the asymptotically constant jump 
across the branch cut of the logarithm. Thus we conclude that the distribution of residues 
should be such that for large n r n should be asymptotically bounded from above. This is 
the condition that the Pade approximant is not satisfying. 

Now, we claim that this condition on the r n in fact prohibits any function where the 
separation between poles grows asymptotically, not just the Pade approximant. For example, 
there is no function with a similar mass pattern m 2 n ~ n 2 but with nonincreasing residues 
that provides a match to PT. This is intuitively clear: adding together functions supported 
on a narrow interval, growing progressively further apart, will never produce a constant. We 
can formalize this: take the average value of pa between two consecutive poles m 2 n and 
If pa is to agree with PT, this average should not go to zero as n — > oo. But we can just 
compute the average: 

- i r l+i i 

PA '" " m 2 n+1 - ml U V ( s - m ^ + A2 ' 
Interchanging the order of summation and integration, this is just 

PA „ = — ^ ? Y ? ( arctan m * ~ f n+l - arctan < ~ < ) . (A.4) 

This expression is valid in general, regardless of the behavior of the m 2 n and r n . 

Now we exploit the properties of the arctangent. If the distance between successive poles 

9 9 

is asymptotically growing, we can choose n such that n+1 A — - ^> 1. Then we use the fact 
that: 

arctan(x) — arctan(a; — Xq) ~ tt6(x)6(xq — x), Xq ^> 1. (A. 5) 

In other words, for very large values of Xo this function is approximately tt on the interval 

9 9 

from to Xo and outside that interval. In particular, in our case we identify x with k A n+1 

9 9 

and Xq with n+ \ — - . It is then clear that the only k values that contribute to the sum in 



Eq. |A.4| are k — n + 1 and k = n. The result is that 



1 (r n + r n+ i) arctan(m2 x - m 2 n ) 
pA,n « — 2 . A.6 

m n+l ~ m l A 

We have argued that the r n are not asymptotically growing; the arctangent is also bounded 
(it is approximately 7r). Hence, if the denominator m^ +1 — m 2 n is growing, this expression 
tends asymptotically to zero at large n, in contradiction with PT. 
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At this point we have established that the radial Regge physics, m 2 n ~ n asymptotically, 
is the behavior which matches PT with the fewest poles on a given interval. We do not have 
an argument that we cannot have slower growth, for instance m\ ~ \Jn. Intuitively it is 
apparent that if we do have slower growth, the residues must decrease or we will construct a 
function that greatly exceeds the perturbative answer at large s. Let's see what happens if 
we attempt to extend the above argument to the case of poles that grow successively closer. 
In that case we are interested in the opposite limit of the difference of arctangents, Xq — > 0. 
This is of course determined by the derivative: 

arctan(x) — arctanfx — Xq) ~ — , Xq <C 1. (A. 7) 

x l + 1 

Assuming the separation between poles is decreasing with n, we can choose n such that 

2 

- <C 1 and the above formula is valid. Note that ^f^- > y provided |x| < 1. Again 



2 n + l 



A 

„2 m 2 



identifying x with k A n+1 and summing only over the terms where |x| < 1, we find that 
the sum is at least 



PA ' nZ ml +1 -ml ^ A 2A = 2A^ ^ ^ (A ' 8) 

n+1 n k: \ml-m 2 k \<A k: |m2-m||<A 

From this we get an estimate of how quickly the residues must decrease: fast enough to 
compensate for the number of poles in a given interval. Assuming ml is a smooth function 
of n, for large n the number of poles in a unit interval is characterized by {ml +1 — m^) _1 , 
so the residue should satisfy r n w ml +1 — m^. 

Now, since we don't have any argument against decreasing poles, it seems possible that 
a function like 



could provide a good match to PT. We have numerically investigated this, and indeed such a 
function appears to agree well. It does not seem possible to determine the correct behavior of 
the distribution of poles just from matching the leading perturbative behavior of two-point 
functions. On the other hand, the asymptotically constant distribution of pole expected from 
Regge physics is in some sense the minimal choice, so in the absence of other arguments it 
seems preferable. 
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B Lorentz Invariance of Flux tube Breaking 

Here we show that for center of mass energies much larger than the the typical transverse 



momentum y/a, the simple model of Section 5.2 gives a Lorentz invariant description of the 



early decays of a flux tube. At large enough center of mass energies the flux tube stretches in 
a straight line between the back-to-back quark and anti-quark, and the transverse momenta 
are negligible. In this limit, the situation reduces to a 1+1 dimensional problem. The 
discussion can be further simplified by noting that the pair production which breaks the flux 
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tube is just like 2D bubble nucleation, which is Lorentz invariant around the center of the 
bubble [60-62]. 

Let us take the string axis to be along the x direction, and our string breaking pair 
production to be centered at x = 0, t = 0, which means that the original production of 
the flux tube took place at some earlier (negative) time. At t = 0, we will take a quark 
of mass m to be at x = xq and moving in the positive x direction with momentum po] the 
anti-quark at the other end of the flux tube is somewhere along the negative x axis moving 



in the negative x direction. The Nambu-Goto-quark equations of motion (5.7) tell us that 
the momentum of the quark will decrease linearly with time: 

Pq =p -at. (B.10) 

When the flux tube breaks, a new quark — anti-quark pair is formed at t — 0. Neglecting 
the small transverse momentum, the new pair is at rest, and in order to conserve energy, 
a section of the flux tube must be removed to account for their rest energy. Thus the new 
anti-quark of mass m is at 

x = d = m/a (B.ll) 
while the new quark is at x = —d, so that the flux tube section from —d to d has donated its 



energy to provide the new masses. Again the Nambu-Goto-quark equations of motion (5.7) 
tell us that the subsequent momentum of the new anti-quark is given by 

Pa = (Tt . (B.12) 

Integrating the velocity, v a = p a / E a , of the new anti-quark we find that its subsequent 
position is 



v a dt= I — Pa dt = Vd 2 + t 2 . (B.13) 



This is, of course, exactly what we expect from 2D bubble nucleation: the trajectory should 
be an hyperboloid, corresponding to a circle in Euclidean space (x 2 + 1\ = d 2 ). Similarily, 
the trajectory of the original quark can be found by integration as well but the result just 
corresponds to another circle in Euclidean space, (x — x c ) 2 + (t E — t c ) 2 = d 2 , however one 
that is not centered about the origin. The center is given by 



a V (J 



r, = .,-„ + v/m2 _ +P ° =x + ^ 2 + 4, (B.14) 



t c = ^. (B.15) 



a 

Thus the position of the quark is simply given by 



x q = x c - ^/d 2 + (t c - t) 2 . (B.16) 
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Now let us calculate the rest mass of the daughter flux tube. This is simple at t = 0. 
The total momentum is just that of the quark (since a string cannot carry longitudinal 
momentum), so p = Pq. The total energy is 



E = m + a(x — d) + W m 2 + p\ = a x + W m 2 + p% . (B. 17) 

Thus the rest mass is 



rarest = \/E 2 —p 2 = ym 2 + o 2 x\ + 2ax ^J m 2 + pi . (B.18) 

What if we had chosen to calculate the rest mass in some other frame? The issue may 
seem quite complex since a boost will change the length of the string, but the string length 
is changing with time so the new time slicing will actually involve the string ends at different 
stages of evolution. However the situation is really quite simple. If we boost to a frame 
moving with velocity v to the right, special relativity tells us that in that frame the energy 
and momentum will be 



E' = >y(E - vp ); p' = r/(po - vE) . 



(B.19) 



Since the new anti-quark trajectory is boost invariant it will be at x' = d at t' = 0, so in 



analogy to Eq. (B.17) we should also have that at t' = the total energy in this frame is 



E' 



ax 



' + ^m 2 + 



P 



/2 



(B.20) 



where x' is the position of the quark in the boosted frame at tf = 0. Comparing Eqs. (B.19) 



and (B.20 ) tells us what x' should be if the model is Lorentz invariant. We can independently 



calculate x' from the boost of the quark trajectory. This is most easily done by rotating the 
center of the corresponding circle in Euclidean space (which is equivalent to boosting the 
the focus of the hyperboloid in Minkowski space). This gives, at t' = 0, 



7 (x c - v t c ) - \Jd 2 + ~i 2 (t c -vx c y 



(B.21) 



which is exactly what is required for Eqs. (B.19) and (B.20) to agree. Thus the rest mass of 



the daughter flux tube is the same in any boosted frame, and the model is indeed Lorentz 
invariant. 
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